Model for Polymorphic Transitions in Bacterial Flagella 
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Many bacteria use rotating helical flagellar filaments to swim. The filaments undergo polymor- 
phic transformations in which the helical pitch and radius change abruptly. These transformations 
arise in response to mechanical loading, changes in solution temperature and ionic strength, and 
point substitutions in the amino acid sequence of the protein subunits that make up the filament. 
To explain polymorphism, we propose a new coarse-grained continuum rod theory based on the 
quaternary structure of the filament. The model has two molecular switches. The first is a double- 
well potential for the extension of a protofilament, which is one of the eleven almost longitudinal 
columns of subunits. Curved filament shapes occur in the model when there is a mismatch strain, 
i.e. when inter-subunit bonds in the inner core of the filament prefer a subunit spacing which is 
intermediate between the two spacings favored by the double- well potential. The second switch is a 
double- well potential for twist, due to lateral interactions between neighboring protofilaments. Co- 
operative interactions between neighboring subunits within a protofilament are necessary to ensure 
the uniqueness of helical ground states. We calculate a phase diagram for filament shapes and the 
response of a filament to external moment and force. 



I. INTRODUCTION 

Allosteric proteins switch between different conforma- 
tions in response to the binding of ligands. Allostery 
underlies a wide range of cellular processes, such as the 
cooperative binding of oxygen to hemoglobin [| , the neg- 
ative feedback control of the DNA-binding trp repressor 
in Escherichia coli |2|, and the regulation of the cata- 
lyst aspartate transcarbamoylase, also in E. coli [3J. All 
these examples involve conformational changes of a small 
molecule in response to the the binding of a ligand. Bray 
and Duke have recently argued that the concept of al- 
lostery should be extended from small molecules to large 
assemblies of proteins, in which conformational change 
propagates through a one- or two-dimensional lattice |^| . 
In this picture of "conformational spread," the state of 
an individual protein subunit of a large assembly is sta- 
bilized either by ligand binding or by the conformational 
state of a neighboring subunit. 

In this article we develop a theory based on confor- 
mational spread for polymorphic transformations of the 
flagellar filaments of bacteria such as E. coli. These fila- 
ments are hollow tubes built from identical protein sub- 
units known as flagellin h\. To explain why the tubes 
are usually helical, instead of straight cylinders, Asakura 
supposed that each flagellin subunit can be in one of two 
distinct conformations, with one slightly longer than the 
other |6j . Since the subunits may be grouped into eleven 
protofilaments which slowly wind around each other to 
form the filament, a helical filament shape arises when 
some of the protofilaments consist of only long subunits, 
and the rest consist of only short subunits. For example, 
for physiological solution conditions and in the absence 
of external forces and moments, the filament of E. coli 
is left-handed with a pitch of 2.5 /jm and a diameter of 
0.4 /im 7]. Indirect evidence of the spread of conforma- 
tional change comes from the response of a helix to hydro- 
dynamic torque, which can trigger a polymorphic tran- 



sition from a left-handed normal state to a right-handed 
state 8]. For example, Hotani has detached filaments 
from Salmonella typhimurium cells and used hydrody- 
namic torque to trigger polymorphic transformations |9( . 
A rotating helix in a viscous fluid generates propulsive 
thrust; conversely, flow directed along the axis of a helix 
not only stretches the helix, but also tends to untwist 
the helix, independent of its handedness. By flowing so- 
lution past filaments stuck at one end to a microscope 
slide, Hotani observed different helical states propagat- 
ing from the stuck end to the free end . 

Polymorphic transitions caused by hydrodynamic 
torque also occur in free-swimming bacteria ||. There 
are usually five to six flagella per cell. When the flagellar 
motors turn counterclockwise (as viewed from outside the 
cell), the left-handed helices wrap into a bundle which 
propels the cell forward 0, |n| . Reversal of the rota- 
tion direction of a motor leads to a sequence of events 
in which the corresponding filament unwinds from the 
bundle and transforms to a right-handed state be- 
fore eventually wrapping up into the bundle once the 
motor returns to counterclockwise rotation. Fluorescent 
labelling of the filaments has recently revealed the pre- 
cise dynamics of this process jjj. The typical sequence 
begins with a transformation from the left-handed nor- 
mal state to a right-handed "semi-coiled" state with half 
the helical pitch of the normal state. This transformation 
changes the swimming direction of the cell. Then there 
is an additional transformation to a right-handed state 
called "curly-1," which has a similar pitch but roughly 
half the radius of the semi-coiled state. When the motor 
reverses back to counterclockwise rotation, the filament 
transforms directly to the normal state. 

Discontinuous transitions among these states can also 
occur due to ch ang es in pH [TJ, salt concentration [l^ . 
or temperature Il3l . Straight states |(j and curved planar 
( "coiled" ) states |T3| have also been observed. 

In this article we describe a continuum model for flag- 
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FIG. 1: The lattices of (a) L-type and (b) R-type filaments, 
at radius a — 4.5 nm and using data from flEf . The vertical 
edges are identified. Note the II-, 6-, and 5-start helices. The 
protofilaments are the Il-start helices. 



ellar filaments. We begin in section [D] with a review of 
experiments that reveal the structure and mechanical re- 
sponse of the flagellar filament, and then discuss previous 
theoretical models for polymorphism. We introduce our 
model in section ITTT1 The essential features of our model 
are two inequivalent subunit conformations, a mismatch 
of the preferred spacing of protein domains in the inner 
and outer core of the filament, and cooperative interac- 
tions between neighboring subunits. All three of these 
features are suggested by the experiments, and are re- 
quired for our model to yield a unique helical ground 
state for given parameters. In section lTVI we study a sim- 
plified version of our model, introduced in which 
neglects extensibility and twist-stretch interactions. Us- 
ing this model, we calculate the filament ground states 
as a function of material parameters. Then we turn to 
the response of the filament to external moment. Some 
of these results were described in 01 • We conclude this 
section with new calculations of the response to exter- 
nal force. In section we present the phase diagram 
and response to external moments and forces for the full 
model, in which the simplifying assumptions of section lTVl 
are removed. We discuss our results in section IVT1 



II. SUMMARY OF PREVIOUS WORK. 

Before the detailed structures of virus capsids were 
measured, it was argued that self-assembly of identical 
subunits would lead to highly symmetric forms, with all 
subunits in equivalent bonding environments 16] . For 
example, the protective coat of tobacco mosaic virus con- 
sists of identical protein subunits which form a straight 
cylinder 0] . In flagella, the flagellin subunits also form 
a two-dimensional crystal on the surface of a cylinder, 
which is about 20 nm in diameter. As mentioned already, 
the subunits can be grouped into eleven protofilaments 
that gradually wind around each other to form the fil- 
ament. Since the flagellum usually has a helical shape, 
exact equivalence cannot apply. However, since the fil- 
ament is thin, small departures from exact equivalence 
can lead to helices of the observed dimensions 61. 




FIG. 2: Cartoon of a short segment of the R-type flagellar 
filament, showing the outer core (A), connecting spokes (B), 
inner core (C), and lateral bonds (D) between neighboring 
protofilaments (E). The position of the spokes is not accurate, 
and is only meant to indicate that the inner and outer cores 
are connected. The outer domains are not shown. The lateral 
bonds lie along the five 5-start helices. 



A. Experiments 

A consequence of Asakura's picture |6| of two stable 
states for protofilaments is that two distinct straight fil- 
aments are possible, one in which all protofilaments are 
short, and one in which all are long. These states, known 
as R-type and L-type, have been observed 18]. The sub- 
unit period along a protofilament is 0.8 A shorter in R- 
type than in L-type. The protofilaments make a positive 
angle with the longitudinal direction in R-type, leading 
to right-handed twist, whereas the protofilaments in L- 
type have a left-handed twist. Figures Q (a) and (b) 
show the subunit lattices in L- and R-type at a radius 
of 4.5 nm, obtained from x-ray scattering from fibers of 
filaments The primitive lattice vectors a and b in 

the two figures are slightly different, but in both cases 
6a + 5b = Cx where C is the circumference of the fila- 
ment at 4.5 nm, and Cx = since the lattices lie on the 
surface of a cylinder. The vector constraint on the two 
lattice vectors yields two degrees of freedom, which may 
be taken to be the subunit spacing along the 11-start di- 
rection and the angle the 11-start direction makes with a 
longitude. Equivalently, we may replace the angle degree 
of freedom with the twist of the filament. 

Electron microscopy has shown that the filament cross- 
section consists of outer domains surrounding a core with 
inner and outer parts (Fig. |2J) lij . The inner and outer 
cores are connected by radial spokes. The outer domains 
do not seem to affect the mechanical properties of the 
filament; for example, mutations in the outer domain do 
not affect polymorphism [20|. Both the outer core and 
the inner core are crucial for polymorphism. Mutations 
in the outer core can change the ground state of the fila- 
ment from normal to curly, semi-coiled, or straight |2l| . 
Without the inner core, there is no polymorphic behav- 
ior: the filaments are straight, and the subunit lattice at 
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radius 4.5 nm is of Lt-type . w ith left-handed twist and 
subunits in the short state [22]. (But note that subunits 
without the inner region can polymerize to form helical 
filaments if there is a wild- type helical seed p3|). 

Since flagellin tends to self-assemble into filaments, it 
is difficult to crystallize. By clipping off some of the 
residues at either end of a flagellin chain, Samatey et 
al. crystallized the F41 fragment of flagellin from the 
R-type filament of S. typhimurium, allowing 2.0 A reso- 
lution of the subunit structure by x-ray scattering |24| . 
The crystals consist of an anti-parallel arrangement of 
R-type protofilaments. The lateral interactions normally 
present in an intact filament arc different in the crys- 
tal, but the protofilament is still stable. Therefore, the 
authors of [24j suggested that the protofilament is an "in- 
dependent, cooperatively switching unit." The structure 
of the intact R-type filament was revealed at the 4 A res- 
olution level by electron cryomicroscopy [25|. This work 
showed that each subunit interacts with its neighbors in 
the outer core along the 5-, f f-, and 16-start directions, 
but not along the 6-start direction. In the inner core, 
each subunit contacts its neighbors along the 5-, 6-, and 
Il-start directions. 



B. Models 

The first quantitative theory of polymorphism in flag- 
ella is due to Calladine, who modeled an element of a 
filament as two rigid discs connected by eleven linear 
springs [2(| H^] . Following Asakura, he supposed that 
the springs could have one of two rest lengths, with the 
same spring constant in either case. If all the springs have 
the same rest length, then a stack of elements would be 
straight. If some are short and the rest are long, then 
in equilibrium there will be a variation of spring length, 
causing the discs in an element to incline relative to each 
other. A stack of these elements will form a curved fil- 
ament in equilibrium, with the curvature varying sinu- 
soidally with the number of long springs [27j • 

Calladine further supposed that the straight state with 
long springs has the twist of the L-state, the straight 
state with short springs has the twist of the R-state, and 
a filament with a mix of short and long springs has a 
twist that linearly interpolates between the two extremes. 
These assumptions lead to ten discrete helical states with 
curvature varying sinusoidally with twist. A key assump- 
tion of this model is that the torsion of a helical state is 
equal to the twist of the elements. Calladine's model was 
studied further by Hasegawa et al. [29j, who used more 
accurate data to determine the twist parameters of the 
model, and who also calculated the range of deviations 
in subunit spacing in single filament of a given helical 
shape. The predictions of Calladine's model for the pos- 
sible ground states of flagellar filaments agree reasonably 
well with measurements of the filament shape [l5| . 

Goldstein et al. used a different approach to model 
flagellar filaments |30| . They presented a continuum the- 



ory for an elastic rod with a preferred direction of curva- 
ture, a preferred magnitude of curvature, and a double- 
well potential for twist. Although the preferred direction 
of curvature is hard to reconcile with the fact that the 
flagellin subunits are identical, this model is probably 
the simplest continuum model for bistable helices. Un- 
like the model of j2(| l2?l |28| , it provides a framework for 
calculating the deformation of the filament in response 
to external forces and moments, such as hydrodynamic 
loading jsj. This continuum model, along with resis- 
tive force theory for hydrodynamic drag, leads to a pre- 
diction for the velocity of the front propagating along a 
helix du ring a transition from one polymorphic form to 
another |31|. 

In this article we develop a new model which has el- 
ements similar to both the spring model and the con- 
tinuum bistable helix model 0]. Like the model of 
Goldstein and collaborators, our model is a continuum 
rod theory that predicts the deformation in response to 
load. Like Calladine's model, our model is based on the 
two conformations of the protein subunits. However, our 
model introduces important new features. For example, 
it does not single out a preferred direction of curvature, 
but instead assigns every protofilament the same double- 
well potential for stretch. The model incorporates pre- 
stress by assuming there is an elastic mismatch between 
the inner and outer cores of the filament; this elastic 
mismatch is necessary to make the absolute minimum 
of the elastic energy helical instead of straight. Finally, 
we explicitly treat the cooperative interactions between 
neighboring subunits. These cooperative interactions are 
required to make the ground state unique. 



III. FILAMENT ELASTIC ENERGY 

The goal of this article is to explain why flagellar fil- 
aments are helical and why they undergo polymorphic 
transformations. Since the characteristic length scale of 
the helical polymorphs is a micron, much larger than the 
5 nm scale of the flagellin subunits, it is natural to formu- 
late a coarse-grained or continuum theory for polymor- 
phism. We model the filament as two concentric elastic 
cylinders (Fig. [2]). The inner cylinder corresponds to the 
inner core of the filament, and is treated as an elastic rod 
with a resistance to stretching. The outer cylinder corre- 
sponds to the outer core of the filament, and consists of 
eleven strands which gently wind around the inner core 
(Fig. [5J. Each strand corresponds to a protofilament, 
although strictly speaking, a protofilament in a flagellar 
filament includes material from both the inner and outer 
cores. We divide the outer core into protofilaments but 
treat the inner core as a single elastic element since only 
the outer core is implicated in switching between poly- 
morphic forms |2l| . 

Our model treats the eleven strands and the inner core 
as continua. To define the strains of the strands and the 
inner core, we choose a reference configuration and label 
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FIG. 3: (a) The two preferred conformations of the subunit. 
(b) Sketch of double-well potential for stretching a protofila- 
ment, with minima corresponding to the two preferred states 
of (a), (c) The two preferred relative positions of neighboring 
subunits. (d) Sketch of double- well potential for twist K3. 

the points along the inner core with arclength S in the 
reference configuration. The strain of the inner core is 
e = ds/dS — 1, where s is the actual arclength of the 
inner core. Note that the strain need not vanish when 
the inner core is in the reference configuration. Since the 
difference in subunit spacing of L- and R-type is about 
2% of the subunit size, we assume e <C 1. To define the 
strain Ci of the ith strand, where i = 1, 2, . . . , 11, consider 
the cross-section of the filament at S. A point on a strand 
that intersects this cross-section is assigned the label S. 
Thus, if d£i is the element of length of the ith strand, 
then a = diJdS- 1. 

Our model treats the inner core as a rod with a 
quadratic strain energy per unit length, 

U = h s (e-e ) 2 , (1) 

where k s is an elastic constant and e is the strain of the 
inner core in its unstressed state. The total stretching 
energy of the inner core is J Uq(S) dS. But since dS = 
ds/(l + e) ~ ds to an excellent approximation, we may 
approximate all energies as integrals over s rather than 
S. The distinction between s and S must be maintained 
only when computing strain e. 

Motivated by the discussion of two preferred subunit 
conformations, we introduce a "switch" for protofila- 
ment extension by defining the stretching energy per unit 
length of a strand as a double-well potential, 

<J SJ = M [( e .?-e2) 2 /4- Ul e^], (2) 

where u is an elastic constant, e p determines the posi- 
tions of the local minima of U S i that correspond to the 
two conformations, and the dimcnsionlcss parameter u± 
determines the degree of asymmetry of U S i [Fig. EJa), 
(b)]. A term proportional to ef is disregarded for sim- 
plicity. 

Since radial spokes connect the inner and outer cores, 
the strains e and are not independent. For example, 
when the filament is straight, e = e;. Mismatch strain 



arises because the preferred strain of the inner core need 
not agree with the either of the strains corresponding to 
the minima of U S i- We will vary the mismatch strain by 
changing eo and keeping U S i fixed. Mismatch strain is a 
natural assumption since during self-assembly, the sub- 
units must denature as they travel up the narrow channel 
at the core of a filament before inserting into their final 
destination at the distal end |l(J . 

When the filament is curved, the extensions e, depend 
on i. To determine this dependence, we follow Calladinc 
and make the simplest assumption by supposing that pla- 
nar cross-sections of the filament remain planar under 
bending and twisting [27]]. We further suppose that the 
rod is unshearable. Thus the position of a material point 
on the ith strand at radius a is 

Yi = r c + a cos /3,-ei + a sin $e 2 , (3) 

where r c is the position of the centerline of the inner 
core, fa = 27r(i — 1) /ll, and ei and e 2 are orthonormal 
material frame vectors which are perpendicular to the 
tangent vector 63 = dr c /ds (Fig. |2J). Without loss of 
generality we have chosen ei to point to the protofilament 
with i = 1. 

In this article we choose a = 4.5 nm, the radius at 
which the lattice structures of Fig. 2] are reported. Since 
{ei,e2,e3} form an orthonormal frame, 

^=KXe„ (4) 

where n — n^e^, n\ and K2 are the components of the 
curvature vector k± = de^/ds = K2<2i — Kitb, and K3 
is the twist. The components K\ and K2 are the rates 
that the orthonormal frame rotates about ed and 62, re- 
spectively. Likewise, the twist k 3 is the rate that the 
orthonormal frame rotates about 63. Note again that 
since eCl, whether we use d/ds or d/dS to define k is 
immaterial. The magnitudes of the components of k are 
much smaller than 1/a: \ny\a <C 1. Thus, the relations 
© and Q together with d£ = (dr^/dS* • dv./dSf^dS 
imply 

d ^ _ dS , -a a «\ 

ti = — — w e + clki sin Pi — aK2 cos/v (5) 

db 

Although the strain energy of Eq. involves terms of 
up to fourth order in e,, the form of the coefficients of the 
quartic potential justifies approximating to first order 
in e and n^a. 

Note that the twist K3 does not appear to first or- 
der since the protofilaments are almost longitudinal. For 
example, at radius 4.5 nm, the angle between a protofil- 
ament and the longitude is —1.45° = — 0.0253 rad in L- 
type and 3.46° = 0.0605 rad in R-type 

The stretching energy of all the strands is U s — 
Sjirf/si, or 
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FIG. 4: Illustration of how twist causes neighboring subunits 
to slide past each other. The subunits rotate rigidly but the 
lattice of centers of mass shears. The angle of shear has been 
exaggerated for illustration. 



where k = \k±\ > is the curvature of the filament. 
Although U s i has a term linear in ej, there are only 
even powers of k in U s . Thus, the energy U s depends 
on the magnitude of the curvature. There is no natural 
(or spontaneous) direction of curvature since all subunits 
are identical. However, the subunits may be in different 
states. For example, if 3e 2 — e 2 < 0, then there will 
be a natural magnitude of curvature, and a distribution 
of subunit states around the circumference of the fila- 
ment. In this case, the direction of the curvature vector 
at a point s may point in any direction in the plane of 
the cross-section. Note also that U s shows no trace of 
the original eleven-fold symmetry of the sum of U S i over 
the eleven protofilaments; the resistance to bending is 
isotropic for a rod with a cross-section with a symmetry 
greater than two- fold [32]. This isotropy rules out cou- 
plings between twist and curvature like those that arise 
in the mechanics of DNA |3^| . 

In our model, lateral interactions between neighboring 
subunits on different protofilaments lead to a resistance 
to twist. We follow the suggestion of Namba and Von- 
derviszt [34| and suppose that neighboring protofilaments 
have two preferred relative displacements along the local 
protofilament direction [Fig. EJc), (d)]. These interac- 
tions are at the quaternary level; they do not require a 
conformational change of the subunit. As Fig. 0] illus- 
trates, twisting a filament causes neighboring protofila- 
ments to slide past each other. Therefore, two preferred 
relative displacements are equivalent to two preferred 
twists. Thus, we introduce a second switch, a double- 
well potential for twist: 



(J) 



rip determines the two 



where v is an elastic constant 
preferred twists, and v\ makes the potential asymmetric. 
As for our potential for extension, Eq. a cubic term 
is disregarded for simplicity. 

In addition to the switches for extension and twist, the 
filament structure suggests a coupling between twisting 
and stretching. As mentioned earlier, the electron cry- 
omicroscopy studies of the intact R-type filament show 
that each subunit has contacts with its neighbors along 



the 5-, 11-, and 16-start directions in the outer core. Our 
model effectively replaces all the inner core lattice con- 
tacts with the inner tube, and those along the 11-start 
direction of the outer tube by a double-well potential. 
Since the 16-start direction is almost longitudinal, like 
the 11-start direction, we ignore it for simplicity. The 5- 
start helices make an angle of approximately —45° with 
the longitude in both the straight filaments. This con- 
tact will lead to twist resistance, as well as a coupling 
between twist and stretch. To understand the sign of 
this coupling, consider a straight R-type filament in the 
limit where the bonds along the 5-start direction are rigid 
and unable to stretch. Pulling on the filament will cause 
the left-handed 5-start helices to untwist a little, lead- 
ing to an increase in K3, the twist of the 11-start helices. 
Likewise, an applied moment which tends to unwind the 
5-start helices will cause the filament to lengthen. Thus, 
the twist stretch coupling is of the form 



U 5 



-k 5 aeK 3 , 



(8) 



where k 5 is a positive elastic constant. 

The 5- and 6-start contacts in the inner core will also 
contribute to the twist-stretch coupling. Without know- 
ing the relative strength of these bonds, we cannot pre- 
dict the sign of the contribution to the twist-stretch cou- 
pling, since the 5-start helices are left-handed and the 6- 
start helices are right-handed. However, since the outer 
core bonds lie at a radius roughly double that of the in- 
ner core bonds, and bending and twisting moduli in rod 
theory scale as the third power of radius for a hollow 
rod [32|, we assume that the contribution to the elastic 
energy of the inner core bonds is dominated by that of the 
outer core bonds. Note also that the 5-start bonds will 
contribute to the bending and twisting resistance of the 
filament; these effects are already accounted for in the 
bending terms of U s and the twist potential U t , Also, 
since the high-resolution structure of the intact L-type 
filament is currently unavailable, we make the simplest 
assumption and suppose there are contacts along the 5- 
start direction but not the 6-start direction in L-type as 
well as R-type. Finally, we note that there could be a 
twist-stretch coupling arising from the elasticity of the 
subunit itself. For example, pulling on a subunit along 
the 11-start direction could cause it to shear as well as 
stretch, causing the whole filament to twist in a right or 
left-handed sense. We know of no evidence for or against 
this possibility, but we must acknowledge that our choice 
of the sign of k§ is provisional. 

The terms of the elastic energy introduced above (t7 s , 
Ut, and U5) determine the extension e, the curvature k, 
and the twist K3. However, these variables do not com- 
pletely determine the shape of the filament: two rods may 
have the same curvature and the same twist, but different 
shapes. For example, Fig.[S]shows circular and helical fil- 
aments that have the same curvature and the same twist, 
and hence the same value of U s + Ut + U5 (assuming e is 
the same for both filaments) . To see what kind of term is 
required to remove this degeneracy, consider the special 
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FIG. 5: Two rods with the same curvature and the same 
twist. The black line traces the path of the head of the vector 
ei. The circular rod is twisted through one turn to make 
K3 = 1/Rc, where R c = 1/k is the radius of the circle. 



wa 
~2 



E 

i=l 



(9) 



where w is an elastic constant with units of energy per 
unit length. The cooperative interaction Eq. © stabi- 
lizes the state of a given subunit, depending on the state 
of the subunits neighbors, in accord with the concept of 
conformational spread . In terms of the filament vari- 



FIG. 6: The basis vectors of the material and Serret-Frenet 
frames in the plane of the rod cross-section. 



ables e, «i, and k%, the cooperative interaction is 



case of a rod bent into a circle. In Fig. [5] the curvature JJ C = 

vector of the circle lies in the plane of the circle, pointing 
toward the center. Since the twist is nonzero, k 3 0, 
the components of the curvature vector in the material 
frame depend on arclength, dKi,2/ds ^ 0. Therefore, 
by Eq. ||SJ|, the extension e, of the ith protofilament also 
depends on arclength. Contrast this example with the 
case of a circular rod with no twist, K3 — 0. In this 
case, K\ = and k 2 = l/R, and the extension ej of each 
protofilament is uniform. The shortest protofilament lies 
a circle of radius R — a, and the longest protofilaments 
each lie on a circle with radius slightly less than R + a 
(since there is an odd number of protofilaments). The 
helix in Fig. [5] is analogous to the circle with no twist; 
the pitch and radius have been carefully chosen to ensure 
that ej is uniform. 

Motivated by these observations, we suppose there is 
an energy penalty when neighboring subunits on the 
same protofilament are in different states. This co- 
operative interaction is implicit in the assumptions of 
Asakura |(j and Calladine [26j, who supposed all sub- 
units in a given protofilament prefer to be in the same 
state. Also, as mentioned in section III Al the integrity 
of the protofilaments in crystallized R-type subunits led 
Samatey and collaborators to suppose there are strong 
cooperative interactions along the 11-start direction |24| . 
In our continuum model, the cooperative energy penal- 
izes nonuniform extension of each protofilament: Thus, 
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( 10 ) 

It is also convenient to write U c in terms of the Serret- 
Frenet basis [23|. Recall that the Serret-Frenet basis of 
a space curve is an orthonormal frame consisting of the 
tangent t = 63 = dr/ds, normal n, and binormal b, 
where 



(11) 



and t is the torsion. (Recall also that unlike r, k\, «2, 
and K3, each of which may have either sign, the curvature 
k is nonnegative.) Since {n, b} and {61,62} span the 
same plane, they are related by a rotation: 




ei = sin fh — cos /b 
62 = cos/n + sin/b. 



(12) 
(13) 



These relations, together with k± = de^/ds, imply 
Ki = — ncosf and k 2 — nsmf (see Fig. [BJ. The rela- 
tion between twist and torsion may be found by writing 
«3 = — ei • de2/ds in terms of the Serret-Frenet basis: 



^3 = df/ds + 



(14) 



U c = 



llwa 
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The cooperative term U c is at its absolute minimum when 
the extension and curvature of the filament are uniform, 
and when the torsion is equal to the twist. Thus, U c 
removes the degeneracy in the filament ground states. 
When r = k 3 [as for the helix of Fig. [3 with radius 
R = k/ (k 2 +k|) and pitch P = 27tk 3 /(k 2 -)-K3)], the angle 
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/ is constant, and therefore the Serret-Frenet frame does 
not rotate relative to the material frame as s changes. 
Without loss of generality we may consider / = ir/2, 
in which case the Serret-Frenet frame and the material 
frame align for all s. Therefore, the protofilament i = 1 
traces out the path r(s) + an(s); this path is the "shortest 
path" on the surface of the helical filament, analogous to 
the path of radius R — a on the surface of a circular tube 
with centerline radius R and tube radius a. 
To summarize, the total elastic energy is 



e./e 1.5 



E 



(U a + U s + U t + U 5 + U c )ds. 



(16) 



The material parameters, such as eo and V\, for example, 
are functions of environmental conditions such as pH and 
temperature. These functions must be determined by a 
microscopic theory and are therefore outside the scope 
of our coarse-grained theory. Therefore, in the following 
sections we will determine the filament shape as a func- 
tion of the material parameters and externally imposed 
moments and forces. 



IV. SIMPLE MODEL: STIFF CENTRAL 
SPRING 

First we consider the simplified model studied in [T3 |. 
in which the stretching stiffness of the inner core is large 
(k s /u — > oo), the nonlinear stretching potential is sym- 
metric (u\ — 0), and the twist-stretch coupling is disre- 
garded (&5 = 0). Since the central core is so stiff, the 
extension of the filament is e = eo, and the total energy 
density (up to an unimportant additive constant) simpli- 
fies to 



U 



U 4 I 2 2\ 2 

-a {k - k q ) 



a 4 , . 



llwa -3 



d/ti 
~d7 



( ^Hl 
\ ds 



(17) 



where u = 33m/8 and n^a 2 = 4(e^ _ 3eg)/3. Note that 
since k s /u — > oo, the effect of making fc 5 nonzero may be 
captured by shifting the value of v\ in the model with 

h = o. 



A. Ground states 

To find the configuration of the filament, we minimize 
the energy J Uds by varying ki, and K3. In this 
subsection, we find the configuration in the absence of 
external moments and forces and as a function of the 
material parameters. In our discussion we will only con- 
sider straight and helical filaments. The cooperative term 




FIG. 7: Protofilament extension ei vs. preferred extension eo 
of the inner core for i = 1 (a), i = 2 and i = 11 (b), i = 3 and 
i = 10 (c), i = 4 and i = 9 (d), i = 5 and i = 8 (e), i = 6 and 
i = 7 (f). The region between the horizontal dashed lines is 
the unstable region. 



has an absolute minimum for constant k\ and K2] there- 
fore, as discussed above we take / = tt/2, or k\ = 0. 
Minimizing U over K2 leads to 



K2 = 

for |e | > e p /\/3 (or Kq < 0), and 

4(4 - 3e^)/3 



2 2 
K 2 a 



(18) 



(19) 



for |e | < e p /V3 (or k% > 0). 

To interpret these results physically, consider a straight 
filament. If | eo | < e p /V3, then the strain a of each of the 
eleven protofilaments lies in the unstable region of the 
double-well potential, where d 2 U s /de 2 < 0. Therefore, 
the straight state is unstable. The filament can lower its 
energy by bending, since in the bent state, the protofil- 
aments on the inside of the curve will have an extension 
close to — e p , and the protofilaments on the outside of the 
curve will have an extension close to e p . Figure shows 

the stretch e, = e + (4/3) 1 / 2 y / £ 2 - 3eg cos[27r(i - 1)/U] 

of the individual protofilaments as a function of eo. Note 
that when the filament is bent, there are always some 
protofilaments with an unstable extension |ef| < V3e p , 
even though the filament as a whole is stable. 

This last observation illustrates the importance of the 
central stretching potential, Eq. QJ. If there were no cen- 
tral core, then k s — 0, as in the models of j2(| and [29]| . 
When fc s = in our model, both e and k 2 are deter- 
mined by minimizing Uq + U s . States with «2 ^ are 
still possible, despite the absence of lattice mismatch, but 
straight states always have lower energy. This metasta- 
bility of curved states when fc s = is consistent with 
the experimental observation that filaments formed from 
subunits with truncated ends are straight, unless the sub- 
units grow on a helical seed filament |23. The cooper- 
ative interactions at the junction between the seed fila- 
ment and the filament consisting of truncated subunits 
can favor curvature. 
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-0.4 0.0 0.4 e n /e„ 0.4 0.0 0.4 » 



FIG. 8: (a) Curvature vs. eo and (b) twist vs. vi for the 
simple model of Section ITVI In (b), the solid lines correspond 
to the states of lowest energy, the dashed lines to metastable 
states, and the dotted line to unstable states. 




FIG. 9: (a) Applying twisting moments to the right-handed 
rod (top) causes it to snap through to a left-handed state 
(bottom), (b) The rod on top has K2 > 0; the rod on the 
bottom has K2 < 0. Applying moments as shown to the top 
rod causes it to snap through. 



The twist is determined by minimizing U over K3, 
which yields 

4 - n 2 pK3 = Vl n 3 p . (20) 

When vf > 4/27, there is only one value of K3 that solves 
Eq. PU)|. This solution is stable since d 2 U/dn 2 > 0. 
When vf < 4/27, there are three solutions. One of these 
solutions has d 2 U /dn\ < 0, and is unstable. Of the two 
stable solutions, the one with the same sign as v\ has 
lower energy. 

Figure EJa) shows the dependence of curvature K2 on 
the preferred stretch eo of the central potential, and 
Fig. E£b) shows the dependence of twist K3 on the asym- 
metry parameter v\ . These figures directly show the con- 
tinuous nature of the transition from straight to bent, 
and the discontinuous nature of the transition from left- 
to right-handed. Note that since eo and v\ are inde- 
pendent parameters, the stretch and the twist are not 
coupled in our simple model. 




FIG. 10: Two helical filaments. The black line traces along 
the path of the protofilament i = 1. Top: a filament with 
K2 < and < 0. Applying axial moments of sufficient 
magnitude in the directions shown causes the filament to snap 
through to a state with k,2 > and K3 > (bottom). 



B. Response to external moment 

The goal of this subsection is to calculate k for a fila- 
ment subject to an external moment M = Mi, where the 
z-axis is along the helical axis for helices and along the 
tangent vector for straight filaments. We begin with a 
qualitative discussion. Consider a straight filament with 
a double-well potential for twist and a single well poten- 
tial for curvature [for example, a filament with k 2 < 0; 
see Eq. <|17[) ]. To simplify the discussion, take the twist 
potential to be symmetric: v\ — 0. Suppose the filament 
is in a state of right-handed twist when M — 0. Apply- 
ing a positive twisting moment will tend to increase the 
twist of the filament. As the positive moment increases, 
the twist of the filament will increase continuously. Ap- 
plying a negative moment to the filament will tend to 
untwist it [Fig. |jja), top]. As the magnitude of the neg- 
ative moment increases, the filament will untwist contin- 
uously until it snaps through to a left-handed twist at 
a critical moment. Further increase in the magnitude of 
the moment will cause the left-handed twist to increase 
continuously. Once the moment is released, the filament 
will remain in a left-handed state [see Fig.[ffia). bottom]. 
This state is the same energy as the right-handed state 
when v~x = (and fcs = 0); when V\ > 0, this state is 
metastable. 

The situation is similar for a filament with a double- 
well potential for curvature and a single well potential for 
twist (for example, a filament with Kg > and £l 2 < 0). 
In the absence of external moment, the shape of this fil- 
ament will be an arc of a circle. Suppose K2 > 0; that 
is, the material frame vector ei points to the center of 
curvature, or equivalently, the protofilament with i = 1 
lies on the inside of the curved surface [Fig. 0Td), top]. 
Applying bending moments with sufficient magnitudes 
causes the filament to snap through to a state of nega- 
tive K2- When the bending moments are released, the 
filament remains in a state of K2 < 0. 

A helical filament with double-well potentials for cur- 
vature and twist has snap-through behavior in both the 
curvature and twist, since a moment M along the helical 
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axis may be resolved into a twisting moment M t along 
the filament tangent vector and a bending moment Mb 
perpendicular to the filament tangent vector. Figure 
shows helical states analogous to those of Fig. E^a) and 
(b). 

We now turn to the quantitative determination of the 
response to an external moment. The variables used 
in the calculation of the phase diagram are intrinsic; that 
is, they depend on the shape but not the orientation of 
the filament. Since an external moment specifies a direc- 
tion in space, it is convenient to use extrinsic variables 
such as the Euler angles, which are defined using a space- 
fixed frame. We will use the convention of Love [36| . 
The material frame at s is generated by first rotating the 
frame {x, y, z} about z through an angle ip to obtain the 
frame e^, then by rotating this frame about e' 2 by 9 to 
obtain the frame e", and finally by rotating the frame e" 
by <f) about e 3 ' to obtain ej [see Eqs. I|A1|) - I|A9I) in Ap- 
pendix^. Our convention is to have the Euler angles lie 
in the ranges 



< ip < 2tt 
-tt/2 < < tt/2 
< (p< 2tt, 



(21) 
(22) 
(23) 



so that the component of e 3 along z is always nonneg- 
ative. The formulas for the material frame vectors in 
terms of Euler angles are given in Eqs. I|A7HA9(1 . To de- 
termine the components of the vector k in terms of Euler 
angles, rewrite Eq. |gj as — e^xde u /ds ■ ex, where 
€ijk is the Levi-Civita symbol, and find 



ki = 9 sin 4> — ip sin 9 cos < 
K2 = 9 cos <fi + ijj sin 9 sin < 
K 3 — (j> + ipcos9, 



(24) 
(25) 
(26) 



where the dot signifies differentiation with respect to s; 
for example, 9 = d9/ds. 

In equilibrium, equal and opposite moments must be 
applied to the ends of the rod. The principle of virtual 
work is therefore 



SE - M ■ [ST(L) - ST(0)] = 0, 



(27) 



where JT^(s) is the angle of the virtual rotation of the 
material frame about e M . For example, the angle of a 
virtual rotation of the material frame about is given 
by the 62 component of the infinitesimal change in ex: 
ST 3 = 62 • Se%. The general relation follows by taking 
cyclic permutations, ^T^ = e^ v \8e u -e\. Using Eqs. (| A7I - 
1911 to calculate 



ST\ = — Sip sin 9 cost 
<5T 2 = Sip sin 9 sin <f> - 
5T 3 = Sip cos 9 + 5cp 

yields M ■ ST — MSip + Mcos 



do sin i 

• cos 6 



(28) 
(29) 
(30) 



Taking the variation of E — J Uds with respect to 
(Sx, S 2 , S 3 ) = (ip, 9, cp) leads to the Euler-Lagrange equa- 
tions 



dU 



d dU d 2 dU 



ds 0Z a ds 2 dZ a 



0, 



(31) 



where a runs from 1 to 3. Since U is independent of ip, 
(p, and (p (see the appendix), these equations simplify to 



dU 
~d9 



d dU d 2 dU 
ds dip ds 2 d-ip 
d_dU_ d 2 dU 
ds ~dd ds 2 89 
_d_dU 
ds 36 



= 



(32) 

(33) 
(34) 



To complete the description of the variational prob- 
lem, we assume natural boundary conditions: <5S a (0) and 
5E a (L) are arbitrary, and <5S a (0) = SE a (L) — 0. Noting 
again that U is independent of (p, we find that these nat- 
ural boundary conditions imply 







8U~ 










s=L 


dU 


d 8U~ 




dip 


ds d^p _ 


s=L 


dU 


d 8U~ 




89 


d 


s 89. 


s=L 






~8U~ 








M. 


s=L 





M 


Mcos( 



(35) 
(36) 
(37) 
(38) 



Exactly the same conditions apply at s = 0. 

The equations H32II38J) are derived without making any 
assumptions on the form of ip, 9, and <j>. We now turn 
to the special cases of straight and helical filaments. A 
helix has a parameterization 



v(s) 



sm( 



1P 



■ sin^s) 



sm f 



cos(^s), s cos 9 



(39) 



with 9 = and ip = 0. For constant 9 and constant ip, 
Eq. H35|) with a = 2 (and the corresponding equation at 
s = 0) imply 8U/89 = (ll/2)wa 3 iP<Psin9 = at s = L 
and s = 0. But 9^0 and ip ^ for a helix; therefore 
(p = at the ends of the filament. On the other hand, 
since U is independent of cp and (p, Eq. I|34() implies 



= va 4 + ip cos 9^j (j) + ip cos 9^j 



' cp + ip cos 9 } -Q 2 



(40) 



is constant. Since ip and 9 are both constant, <f> must be 
constant. Therefore, <p = for all s. Note that without 
the cooperative term, this conclusion would not follow. 
With loss of generality we choose <f> = ir/2. A shift in <p 
amounts to redefining which protofilament has i = 1. 



10 



In general, the Euler angle <p differs from the angle 
/ defined in Fig. However, there is a simple relation 
between (p and / for a helical rod. To derive this relation, 
use Eqs. l|A4HA6jl to show 



(tii = deg/ds = 6e" + tp sin 9e 2 



(41) 



Thus, for a helix for which (9 = 0, the normal is parallel 
or antiparallel to e 2 ; n = e 2 and f — <f> when k 2 > 0, 
and n = —e 2 and / = (p + it when k 2 < 0. With our 
choice for <f>, / = ir/2 when n 2 > and / = 3ir/2 when 
k 2 < 0. 

Returning to the Euler-Lagrange equations, we note 
that Eqs. $Q and (EHJ imply dU/d<p = M cos 6» for all 
s. Likewise, since dU/dtp = when 9 = ^ = 0, Eqs. 1-iL'l) 
and (EHJl imply <9l7ch/> = M for all s. The other Euler- 
Lagrange equation, Eq. (KS.''>[1 with = tp = <p = 0, does 
not yield an independent equation. Thus, 



xp sin 6* (tp 2 sin 2 # — Kn | = - — - si i i f 
^cos6»^ 2 cos 2 6»-^) = 



M 
M 



■ cos 6*. 



(42) 
(43) 



Note that under the assumption of uniform n\ and k 2 , 
Eqs. (|42II43|) also follow from the constitutive relations 



M„ 



dU 



(44) 



with Mi = 0, M 2 = M s'md, and M 3 = M cos 6>. 

First we analyze a rod that is helical when M = 0; 
hence k 2 , > 0, and we may assume 9^0 and 9 ^ ir/2. 
It is convenient to define a 2 = Kq + S! 2 , the pitch angle 
#o for zero applied moment, tan#o = Ko/f2 p , and the 
dimensionless moment m — M(u + v) /(uva 4 a 3 ). The 
angle 9 may be eliminated from Eqs. I|42|) and (I43|l to 
yield a simple cubic equation for tp: 



ip 3 /a 3 -ip/a 



(45) 



The curve relating ip/a to m is precisely the same as the 
curve in Fig[8{b). Equation 1451 determines the handed- 
ness of the filament since K3 = ipcos9, and cos# > 0. 
Combining Eq. igSJ) with Eq. $¥2$ leads to 



sin u = 



sin Uq — — 



The sign of k 2 = tp sin 9 is not determined by these equa- 
tions since Eq. I|46|l is invariant under 9^—9. The con- 
sequence of this invariance is that for every M, there are 
two helical solutions with the same curvature and twist, 
one with k 2 < (as in Fig. ^2 top), and one with the 
k 2 > (as in Fig. 1111 bottom). This symmetry remains 
valid in the full model also. 

Under the restrictions of uniform ip, uniform 9, and 
<p> = it/2, these solutions are stable if and only if the 
Hessian matrix 



H 



d 2 U/dtp 2 d 2 U/(dtpd9) 
d 2 U/(d9dip) d 2 U/d9 2 



(47) 




FIG. 11: Two helices with the same curvature and twist, one 
with K2 < (top), and one with «2 > (bottom). In each 
case, the black line traces the path of the protofilament i = 1. 




-1 1 m 



FIG. 12: (a) Signed curvature H2 vs. m for u = v and 80 = 
7r/6 (blue curves) and Go = n/3 (red curves). The dashed 
lines correspond to left-handed states, (b) Twist K3 vs. m 
with parameters as in (a). 



is positive definite, or equivalently, dip/dM > 0. To see 
why these conditions are equivalent, note under the re- 
strictions just mentioned the Euler-Lagrange equations 
are 



dU/dip-M = 
dU/d0 = 0. 



(48) 
(49) 



Differentiating these equations with respect to M and 
solving leads to 



dtp 
dM 



-777 =ff 22 /detH 



(50) 



The matrix element H 22 — 2(u + v)a tp 4 sin 2 9 cos 2 8 > 
as long as 9 7^ and 9 7^ tt/2. Therefore, the solution 
is stable when dtp/dM > 0, or tp 2 > a 2 /3, and unstable 
when dtp/dM < 0, or tp < cr 2 /3. An alternate route 
to this conclusion is to use Eq. H4bfl to eliminate 9 from 
the total potential energy per length U = U — Mtp. The 
energy U in terms of tp is a simple double-well potential 
with bias given by m: 



U 



,4 4 



III! 



U + V 




tP 



(51) 
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We now examine the behavior of the solutions. Sup- 
pose the helix is left-handed with k 2 > for m = [Pi in 
Fig. I12f a) and (b)]. Eq. 146(1 has no solution for 9 when 
ip 2 < ip 2 , where 



±o l 1 



(52) 



for sin 9q ^ 



«s v/(u + v). When sin 2 O < v/(u + v), 9 
decreases to a minimum value of as ip decreases to ip c ' } 
otherwise, 9 increases to a maximum value of 7r/2 as ip 
decreases to ip c . If > f 2 /3, all these solutions are 
stable. This condition places a further constraint on 9q, 
leading to two regimes. In the regime with < sin 2 9 < 
(2/3)v/(u+v), an applied moment to > causes the helix 
to deform continuously into a straight twisted filament 
at to = m c (P 2 in Fig. 112(1 : as the moment increases 
beyond m c , there is a discontinuous transition to a right- 
handed helical filament (P3 in Fig ll2(l . In the regime with 
(v+u/3) / (u+v) < sin 2 9$ < 1, an applied moment m > 
causes the helix to deform continuously into an open coil 
with no twist at to = to c (Q2 in Fig. 112(1 : as the moment 
increases beyond to c , there is a discontinuous transition 
to a right-handed helical filament (Q3 in Fig. 1121) . The 
critical dimcnsionlcss moment is 



U + V . 2 

to c = sin 9q 4 



'± 1 



(53) 



(again for sin 2 9q ^ v/(u + v)). 

If ip 2 < cr 2 /3, or equivalently, (2/3)v/(u + v) < 
sin 2 9o < {v + u/3)/(u + v), then as to increases, the he- 
lix becomes unstable before either the straight or coiled 
state is reached. Thus, there is a discontinuous transition 
directly from a left-handed helix to a right-handed helix 
at to = m c = 2/(3 v / 3) (Fig 

The solution curves of Figs. 1121 and 1131 are symmetric 
under the combined operation m i— ► —to and K3 » — K3, 
since the nonlinear twist potential is even in K3 in the 
simple model with v\ = 0. The full model has v\ 7^ 
and therefore will not have this symmetry. 

Now suppose the filament is straight when M = 0. 
In this case Kq < 0, and there is a branch of solutions 
corresponding to straight filaments with 9 = and ip 
given by 



4> 3 /cr 3 - W>/cr) cos 2 9 Q = mu/(u + v) 



(54) 



The Hessian matrix for the straight solutions is diago- 
nal, with 



H n = m 4 (3^ 2 -cr 2 cos 2 6»o) 

tt i (2/~ 2 • 2 n 1 2 2 

"22 = a ip \uo sin t/g + va cos i 



f0 



(55) 
(56) 



Thus, the straight solutions are stable for ip < iph, where 



1P- 



■ sin Uq 



(57) 




FIG. 13: (a) Signed curvature «2 vs. m for u = v and 8q = 
43° ~ 0.23 897T . The dashed lines correspond to left-handed 
states, (b) Twist K3 vs. m with parameters as in (a). 



Equivalently, the straight states are stable for to < TOh, 
where 



U + V O/i 

TOh = sin 6*0 \; 1 



-~^sin 2 9 . (58) 



The green curves of Fig. ll4l show K2 and K3 as a function 

of TO. 

Once |to| is sufficiently large, there is also a branch of 
helical solutions, with ip given by 



ip 3 /a 3 - cos(29 )ip/<J = to, 



and 9 given by 



v 



sin f = 



v cos 2 9q + u sin 2 9q <t 2 



u + v 



u + v 



1P 2 



(59) 



(60) 



The same stability analysis described earlier, Eqs. 1(481 - 
150(1 , implies that these helical solutions are stable if and 
only if dVVdAf > 0, or i> 2 > (er 2 /3) cos(2<9 ). Using 
Eq. ((60(1 to eliminate 9 from the energy leads to the po- 
tential per unit length 



U 



4 4 
a a 



u sin 4 9, 




uv ip 

— TO — 

U + V <J 



(61) 

To determine the critical m at which helical solutions 
are possible, note that Eq. 1(601) has no solution for 9 
when ip 2 < ip 2 . Therefore, the branch of helical solu- 
tions begins when tp = ip\ and to = TOh, precisely where 
the straight states become unstable. The red curves of 
Fig. 1141 show K2 and K3 for u — v and 9q = tt/6. Note 
that since ip 2 > (c 2 /3) cos(26>o), the branch of helical 
solutions is stable. 

Care must be taken in the interpretation the graphs of 
Fig. [21 A straight filament under external moment can 
destabilize via the familiar twist instability of an ordinary 
elastic rod 32], which is different from the polymorphic 
transitions studied in this article. For a rod with ends 
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FIG. 14: (a) Signed curvature ki vs. m for u = v and So = 
7r/6. The dashed lines correspond to left-handed states. The 
red curves correspond to helical solutions; the green curves 
correspond to straight solutions, (b) Twist k$ vs. m with 
parameters as in (a). 



that are not clamped, the critical moment for this insta- 
bility is Mt w i s t = 2ttA/L, where A is the bending stiffness 
and L is the length. For our rod, A is the coefficient of k| 
in the energy; that is, A = a 4 a 2 u sin 2 9q. Thus, we may 
estimate the length Lh at which a straight rod subject 
the moment Wh just sufficient to allow helical solutions 
will undergo the twisting instability: 



L h = 2 



UK 2 ). 



(62) 



Since kq and f2 p are of the order of microns, L\ Y is about 
6 microns. When L > L^, we expect a straight flagel- 
lar filament will buckle when subject to a dimensionless 
moment less than mt. 



C. Response to external force 

With the simplified limit of a stiff central spring still 
valid, we study the response of the filament to an applied 
external force Fz. As mentioned earlier, we only consider 
filament shapes which are either straight or helical. The 
line of action of the force coincides with the axis of the 
helix, as if there were rigid bars perpendicular to z at 
each end of the helix (Fig. H31l. Thus, the applied force 
also leads to an applied moment. If r(0) is the origin for 
the moment, then the applied moments at the ends of the 
rod are r(0) x Fz = -RFtp(0) at s = and RF<p(0) at 
s = L, where ip (not to be confused with the Euler angle 
cj>) is the azimuthal coordinate and R is the radius of the 
helix. To see that this loading condition is compatible 
with helical shapes, consider the equilibrium equations 



M + r x F = 
F = 0, 



(63) 
(64) 



where M(s) is the moment due to internal stresses acting 
on the rod through the cross-section at s, and F(s) is the 



c 



M=RF r ^) 




FIG. 15: The line of action of the force coincides with the 
axis of the helix, leading to moments applied at each end. 



force acting on the rod through the cross-section at s [22| . 
Since F = Fz is constant, we can integrate the moment 
balance equation 1631 to find M = — r(s) x F + C, where 
C is an integration constant. Imposing the boundary 
condition M(0) = -r(0) x F = RF<p(0) yields M = 
RFip(s) for a helix of radius R. Since the moment has a 
constant magnitude, the constitutive relations (|44|) imply 
that helical solutions with constant curvature and torsion 
are possible. If the line of action of the external force 
passes through the end of the rod rather than the z axis, 
then |M| is not uniform, and the curvature and twist are 
not constant. 

The principle of virtual work for our loading condition 



5E - F • [Sr(L) - 5r(0)] 

- [M(L) ■ ST{L) - M(0) • <5T (0)] = 0, (65) 

where M(s) = RFip(s) and ST^ = (,_ LU \8e u -e\ as defined 
before. To express Eq. (|65|l in terms of the variables ip, 
9, and </>, first use z • [r(s) — r(0)] = L e 3 • zds to write 



F • [Sr(L) - <5r(0)] = F [ sin 659ds. 

Jo 



(66) 



Then use the relation cp — — cos 9&2 + sin 9e-s and 
Eqs. and to write M • ST in terms of ip, 9, and 
(j>. Just as in the case of an applied external moment, the 
Euler-Lagrange equations eventually imply </> = 0, and 
we take <j> = ir/2. Thus, 



[M(L) • ST(L) - M(0) ■ <5T(0)] 



F 



sin 2 9 



(67) 



Carrying out steps analogous to those that led to 
Eqs. H32I - I38J) and specializing to the case of a helix leads 
to 



Ip COS I 



ip 2 (i> 2 sir 



F 



■ cos 9 



2 cos 2 I 



0„ = 



F 



■ sin 



(68) 
(69) 



Note that these equations also follow directly from the 
constitutive relations of Eq. I|44fl. M\ — 0, M 2 — 
-RF cos 9, and M 3 = RF sin 9. 

We were unable to find an analytic solution to Eqs. I|68|) 
and H69|) . Instead, we used a numerical continuation 
method to generate the solutions ip and 9 as a function 
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FIG. 16: (a) Signed curvature K2 vs. F for u = v, and 9o = 
7r/6 (blue curves) and 9q = 7t/4 (red curves). The dashed 
lines correspond to left-handed states, (b) Twist kz vs. F 
with parameters as in (a). 



of F [37J . The advantage of the continuation method is 
that it can negotiate turning points such as those in the 
graph of K3 vs. v\ in Fig. El The results of our numerical 
calculation are shown in Fig. 1161 where we have intro- 
duced the dimensionless force F = F/(ua 4 a 4 ). FigurclTol 
shows that uniform helical solutions only exist in a finite 
range of F. To see that Eqs. (168(1 and (|69[) do not have 
solutions for large F, consider their asymptotic form at 
large F: 



ip' ~ — t-jCOS( 
ua 
p 

ib cos 9 ~ — r sin 2 9. 
wa 4 



(70) 
(71) 



The ratio of these two equations simplifies to tan 2 9 ~ 
—v/u, which is impossible since u > and v > 0. We 
speculate that the filament takes a non-uniform shape 
with varying curvature and twist for values of F outside 
the allowed range of Fig. ^| Note that for a given force 
F and a given signed curvature k 2 , there are both left- 
handed and right-handed filaments of same magnitude 
of twist owing to symmetry of the energy function in the 
simplified version of the model; therefore the solid and 
dashed lines overlap in Fig. I16f a). 

We do not discuss the case of an otherwise straight 
filament subject to external force, since in this case we 
found no examples of transitions from straight to helical 
states. 



V. FULL MODEL 

In this section, we consider the full model, which ac- 
counts for the inextensibility of the filament and the lat- 
eral bonds between protofilaments. We follow the same 
analysis used for the simple model. The twist-stretch 
coupling arising from lateral bonds, Eq. ©, makes the 
phase diagram for ground states in the full model qualita- 
tively different from the phase diagram for ground states 



in the simple model (Fig. 2 of 14]). However, the re- 
sponse of the filament to external moment and external 
force in the full model is very similar to the response in 
the simple model. 



A. Ground states 

The full model includes all the terms of the energy ll(jfl , 
which we write as E — J Vds, where V = Uo + U s + Ut + 
U5 + U c and we note again that to a good approximation, 
dS w ds. As in the analysis of the simple model, the 
cooperative term is minimized when k± = 0. Assuming 
helical or straight states and minimizing V over e, k 2 , 
and K3, we find 



Ci£ 



- c 2 
3e 2 K 2 



&57«3 , 37 2 e«;| 



fin 



_K2_ 

2fi n 



2e p fi2 
8fi p 



llu 



11 

fin 



(72) 
(73) 
0, (74) 



where ci = 2k s — 1, k s = £; s /(22we 2 ), c 2 = 2fc s e /e p + Ui, 
mi = ui/(ue p ), k 5 = /c 5 /(llwe 2 ), and 7 = afi p /e p . In 
this section, we choose representative values 7 = 1.8, 
u = v, fcs = 1 and u\ = —0.5 to illustrate the pre- 
dictions of our model. When k 5 ^ 0, there will be 
no discontinuous transitions between helices of opposite 
handedness if k s is too small; hence, we choose the rel- 
atively large value k s = 8. The two parameters v\ and 
C2 are the molecular switches of the full model, just as 
V\ and eo/e p are the switches of the simple model: pos- 
itive v\ favors right-handed twist, and positive c 2 favors 
longer protofilaments. Equation l|73(l shows that there 
are straight solutions with K2 = and curved solutions 
with K2 7^ 0, just as in the simple model. Evaluation of 
the Hessian using the independent variables e, K3, and k 2 
for the straight states shows that they are unstable for 
e 2 < e p /3. Solving Eqs. (JZ2J) and JJU) for vi and c 2 when 
k 2 = and e = ±e p /\/3 yields the boundaries between 
the straight and the curved states shown in the phase 
diagram of Fig. El The phase diagram also shows the 
regions of metastability (the purple, blue, yellow, and 
brown regions), where there is more than one equilib- 
rium solution to Eqs. ifT^HT^jl . We found the boundaries 
of these regions by brute force; we systematically var- 
ied Vi and c 2 , solving the equilibrium equations for each 
pair of values and noting when the number of solutions 
changed. The twist-stretch coupling k^ tilts the total re- 
gion of metastability, whereas in the simple model the 
region of metastability is vertical, v\ < 4/27 (see Fig 2 
of [Ti|])- The phase diagram also has a dotted line show- 
ing where the two possible solutions of left-handed and 
right-handed filaments have equal energy. 

The full model yields many possibilities for transitions, 
including discontinuous transitions between straight 
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FIG. 17: Phase diagram for 7 = 1.8, u = v, k s = 8, ks — 1 
and Mi = —0.5; vi and C2 are dimensionless. The inclined 
dotted line defines the boundary where the energy of the left- 
handed state equals the energy of the right-handed state. The 
diagonal dash-dot line traces out the values of vi and C2 used 
to calculate the curvature vs. twist plot of Fig. 1181 The points 
X and Y in the metastable regions represent the corresponding 
parameters used for evaluating the response of a helical and 
straight filament to external loading respectively. 



states of opposite handedness, discontinuous transitions 
between helical states of opposite handedness, and con- 
tinuous transitions between straight states and helical 
states of the same handedness. These transitions are 
present in the simple model. The full model also pre- 
dicts discontinuous transitions from a straight state to a 
helical state of opposite handedness, which do not occur 
in the simple model. 

As mentioned above, changes in solvent condition 
change the values of the material parameters. To illus- 
trate this effect, in Fig. ^| we plot k vs. K3 for values of 
(v\,C2) along the dashed line of Fig. El There is a si- 
multaneous jump in curvature and twist of the filament. 
Note that k rises so steeply with twist that k may appear 
to be undergoing an abrupt jump with small changes in 
parameters. We also point out that the value of 7 deter- 
mines the maximum pitch angle that can be attained by 
the helical filament. Changing the value of 7 = afi p /e p 
amounts to changing the relative positions of the min- 
ima in the nonlinear potentials @ and J7J . If 7 is small, 
there is a wider range of possible pitch angles. We chose 
7 = 1.8 for Figs.ll7land ll8l which gives a maximum pitch 
angle of about 35°. Higher pitch angle polymorphs such 
as curly I and curly II require lower values of 7. 



B. Response to moment and force 

We now consider the response of a filament to external 
moment and force in the full model. Since the filament 
is extensible, there is an additional contribution to the 
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FIG. 18: Curvature vs. twist, in dimensionless units, for 
7 = 1.8, u/v = 1, k s = 8, ks = 1 and u\ = —0.5. The 
parameters vi and C2 vary along the dash-dot line of Fig. 1171 
The dotted lines corresponds to metastable states for which 
the lower energy states are of opposite handedness. 



expression (|66[) for the work done by an external force: 

F ■ [Sr(L) — Sr(0)] — —F f (sin 059 - cos 6Se) ds, (75) 
Jo 

where we have used the approximation dS » ds. In addi- 
tion to the Euler-Lagrange equations derived for external 
moment in the simple model, Eqs. I|32H38[1 . suitably mod- 
ified to include both an external moment and force, there 
is a new equation corresponding to the stretch degree of 
freedom e, 



dv ddv 

— - — — - Fcosfl = 0, 
oe ds oe 



as well as natural boundary conditions 



~dV~ 




~dV~ 




.9e_ 


s=0 


<9e_ 


s=L 



= 0. 



(76) 



(77) 



The argument leading to tf> = in the simple model ap- 
plies here as well, leading to 



e e 

-o + ci c 2 - 

ipsrnO ( fc 5 7e 



fc57 , cos 9 



91 



fin 



52 cos ( 



92 + 93 cos ( 



Fe p cos9 (78) 
-F sin 9 (79) 



cos 9 + g 3 sin 2 9 = M 



where F = F/(llue*), M = Mfi p /(llite£), and 
g x = (37 2 e^ 2 sin 2 ^)/(2e p fi2) 



.92 
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ijj cos 9 I ip 2 cos 2 
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FIG. 19: (a) Signed curvature K2 vs. M for u = v, k s = 8, 
&5 = 1, Ca = 4, Si = —0.25, and wi = —0.5. The dashed lines 
correspond to left-handed states, (b) Twist ks vs. M with 
parameters as in (a). 



FIG. 20: (a) Signed curvature K2 vs. M for u = v, k s — 8, 
fcs = 1) C2 = 15, Ci = —0.25, and u\ — —0.5. The red and 
green curves indicate the helical and straight filaments, re- 
spectively. The dashed lines correspond to left-handed states, 
(b) Twist «3 vs. M with parameters as in (a). 



We use the numerical continuation method to solve these 
equations for the cases of applied moment with F = 0, 
and applied force with M = 0. We use the values of 
V% and C2 corresponding to the point X in Fig. 1171 and 
the same parameter values for 7, u, v, k s , fcs, and u\ 
used in section IIVI The point X corresponds to a sta- 
ble left-handed normal helix, and a metastable right- 
handed helix; we study the response of the left-handed 
helix. Figure El shows that the left-handed helical fil- 
ament changes its shape smoothly with increasing mo- 
ment, going through helical metastable states until is 
straightens out, and then abruptly changes to a right- 
handed filament. Note the qualitative similarity with the 
predictions of the simple model, Fig. 1121 An important 
difference is that the asymmetry of the potential in the 
full model favors one handedness over the other; for ex- 
ample, at M = 0, the right-handed solution of the full 
model has higher energy than the left-handed solution. 

The response to moment of a filament which is straight 
when M = is shown in Fig. 1201 The parameters V\ and 
C2 are chosen to correspond to the point Y in Fig. El 
which represents a stable right-handed straight filament. 
At these parameter values there is also a metastable left- 
handed straight filament. As the moment increases from 
zero, the twist increases until eventually the straight fil- 
ament destabilizes at a critical torque, and there is a 
continuous transition to a right-handed helix. When the 
moment decreases from zero, the straight filament un- 
twists and and a left-handed helix forms at a critical 
torque. Due to the asymmetry of the potential for these 
parameter values, the interval of moments for which left- 
handed straight filaments exist is a subset of the inter- 
val of moments for which right-handed straight filaments 
exist. Nevertheless, Fig. [201 is qualitatively similar to 
Fig.IH 

Finally, Fig. shows the response of the left-handed 
filament at X to an applied force, showing the same qual- 
itative features we found in the simple model, Fig. 1161 
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FIG. 21: (a) Signed curvature K2 vs. F for u = v, k s = 8, 
ks = 1, C2 = 15, vi = —0.25, and ui — —0.5. The dashed 
lines correspond to left-handed states, (b) Twist K3 vs. F 
with parameters as in (a). 



VI. DISCUSSION AND CONCLUSION 

We have introduced a new coarse-grained model for 
polymorphic transitions in bacterial flagclla. The key 
elements of the model are a double-well potential for 
protofilament stretch, a double-well potential for twist 
arising from lateral interactions between neighboring 
protofilaments, elastic mismatch between the protofila- 
ments and the inner core, and cooperative interactions 
between subunits on the same protofilament. Motivated 
by the chiral structure of flagellar filaments, we have also 
included a potential energy for twist-stretch coupling in 
the full version of our model. This potential energy is 
in marked contrast with the tight coupling between fil- 
ament twist and protofilament stretch assumed in j2(| 
and [i^l, which leads to a discrete set of possible values 
for the curvature and twist of a filament in the absence 
of loading. Our softer coupling in the full model (fc 5 7^ 0) 
and complete decoupling (feg = 0) in the simple model 
lead to a continuous variation of curvature and twist as a 
function of the material parameters, with discontinuous 
transitions only at special critical values. Our prediction 
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calls for new single-molecule experiments in which the 
shape of a single filament is observed as solvent condi- 
tions are carefully varied. 

The main result of our work is the calculation of the re- 
sponse of a polymorphic filament to moments and forces 
applied at the ends. In the simplified version of our the- 
ory, the assumption of a stiff central core leads to semi- 
analytic formulas for the shape as a function of moment. 
These formulas allowed us to reveal a host of possible 
distinct scenarios, such as the moment-driven continuous 
unwinding of a helical filament to a straight shape, fol- 
lowed by an abrupt transition to a helix of opposite hand- 
edness as the moment increases (Fig. blue curves). 
Another important generic prediction of our model is 
that the curvature and the twist both jump at the same 
critical value of moment. Although our model makes 
clear predictions for the qualitative behavior of filaments, 
a critical requirement for testing the quantitative valid- 
ity of our theory is the determination of the values of the 
material parameters such as u, v, v±, and eo. These may 
be measured by bending and twisting straight filaments, 
or perha ps m ay be computing using the techniques devel- 
oped in [38j. We can make a rough order-of-magnitude 
estimate of the moment required to cause a polymorphic 
transformation by observing that the critical moment in 
the simple model is approximately M u<i 4 Kq, where we 
have assumed for simplicity that u w v, kq m fl p , a w Kq, 
and we have also dropped overall numerical factors (such 
as the distinction between u and u). Since the bending 
stiffness A « u<2 4 Kq, we may write M m Akq for the char- 
acteristic moment. Using ko « 1 pm, A w 10 _24 Nm 2 
(see [3|J and references therein) yields a characteristic 
torque for polymorphic transition of 10~ 18 Nm. This 
torque is comparable to the hydrodynamic torque expe- 
rienced by a rotating flagellum, which we can estimate 
by computing the the torque on the cell body. Assuming 
the cell body rotates at u = 10 Hz ^(j and is a sphere of 
radius R = 1 /im leads to a hydrodynamic torque of the 
order of M = 8mjR 3 Lu w lCr 18 Nm 41]. These values 
match well with those of Hotani, who estimated the or- 
ders of magnitude for the transformation of a normal to 
semi-coiled filament to be 10 -18 Nm, and 10 _19 Nm for 
the reverse transformation |9j . 

The shape-moment curves predicted by our theory may 
be readily tested by new micromanipulation experiments 
in which a filament is subject to prescribed moments or 
forces. It would also be interesting to vary the solvent 
conditions while doing these measurements. Also, since 
our theory makes predictions about the sign of k%, it 
would be useful to experimentally track the sign of «2 
during a transition. For example, for a transition such 
as is depicted in Fig. ICT a). a bead which is stuck to a 
protofilament which is on the inside of the helix before 
the transition would end up on the outside afterwards. 

Our theory also suggests new directions for theoretical 
work. In this article we have considered boundary condi- 
tions of fixed moment or force. If instead we fixed the po- 
sitions and orientations of the ends of the filament ( "hard 



boundary conditions" ) , then the filament could take on a 
shape in which two different polymorphic states coexist 
(see Fig. 1(a) of for a similar experimental example 
with free filaments with subunits of mixed type). We 
could use our theory to analyze these states, generalizing 
the calculations for planar beams with a bistable poten- 
tial for curvature carried out in |42| . Another important 
generalization of our theory would be to study the dy- 
namics of polymorphic transitions, which has been stud- 
ied in the context of a different continuum model in |30| 
and [3l]]. Finally, it is natural to attempt to go beyond 
our protofilament-based picture and treat every subunit 
on equal footing, and try to adapt ideas developed for un- 
derstanding martensitic transformations to polymorphic 
flagella Hf. 
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APPENDIX A: EULER ANGLES 

The sequence of rotations defining the Euler angles is 
described in section llV Bl Carrying out these rotations 
on the initial frame {x, y, z} leads to 

= cos tpx + sin tpy (Al) 

e' 2 = — sin tpx + cos tpy ( A2) 

e 3 = z, (A3) 

then 

e" = cos tp cos 9x + sin tp cos 8y — sin 8z ( A4) 

e'2 = — sin tpx + cos tpy (A5) 

63 = cos ■i/' sin #x + smtpsmOy + cos6>z, (A6) 

and finally 

ei = (— sin tp sin <f> + cos tp cos <f> cos 8)x 

+ (cos tp sin <j> + sin tp cos <p cos 9)y 

~ cos (p sin 8z (A7) 

62 = — (sin ip cos <p + cos tp sin cp cos 6>)x 

+ (cos tp cos cp — sin tp sin <p cos 9)y 

+ sin^sin^z (A8) 

£3 = cos -0 sin 6>x + sin ?/; sin 6>y + cos6>z. (A9) 
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2 sin i 



Plugging Eqs. (|A7HA9(1 into Eq. Q, de M /ds = n x e M , + 2 sin 

yields Eqs. (|24H2bj) . which show that the k m depend on 
9, <f>, ip, and <p. Furthermore, 

K 2 = 6 2 + sin 2 61P 2 (A10) 



<6 U - cos 6ip\ $ + sin 2 6»^ 2 (All) 



and Therefore, k and both independent of ip, 



= sin 2 ^ 2 



2 and (p. Since K3 is also independent of these variables 

ip 2 + 8 1 (j) — cos 6*i/) j [see Eq. I|26|l ]. the energy density U is also. 
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